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In the last ten years, the generalized central limit theorem established by Paul Levy 
in the thirties has been found more and more relevant in physics. Physicists call 
'Levy flights' random walks for which the probability density of the jump lenghts 
x decays as l/x 1+a with a < 2 for large x. We give here a glimpse of Levy 
flights in physics through two examples, without going into technical details. We 
first introduce a simple toy model, the Arrhenius cascade. We then present an 
important physical process, subrecoil laser cooling of atomic gases, in which Levy 
flights play an essential role. 
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1 Introduction 



The 'usual' central limit theorem (CLT) is an essential tool in physics and in other sciences. 
Indeed, one often knows the probability density f{x) of a quantity x associated with single 
events. Then, one wants to derive the probability density /(Xjv) for a sum 



N 
i=l 



(1) 



of a large number N of such events, considered as independent]]. In simple terms, the 
usual CLT tells that the density of (Xjy — N(x))/y/~N tends to a gaussian distribution at 
large N and that this distribution is determined only by the average value (x) and the 
second moment (x 2 ), provided that these quantities are finite^. The condition for applying 
the usual CLT (finiteness of (x 2 )) is so frequently satisfied that most physicists implicitely 
believe that this theorem applies universally. 

However, physical phenomena can exhibit statistical properties that are beyond the 
usual CLT. In particular, densities f(x) with power law tails: 



for 



oo 



(2) 



(with a > to ensure normalizability) are simple laws that tend to appear frequently. 
If a > 2, (x 2 ) is finite and the usual CLT applies. On the contrary, if a < 2, (x 2 ) 
diverges^] and the usual CLT does not apply. If a < 1, even (x) diverges. As stated by 
the 'generalized' CLT, if < a < 2, the density of Xn (properly renormalized) still tends 
to a stable law, which is not a gaussian but a Levy law. After Mandelbrot, we call 'Levy 
flight' a random walk in which the probability density of jump lengths is given by Eq. @ 
with a < 2. 

The generalized CLT was already available in the thirties but, surprisingly, it has had 
a limited influence on physics for a long time. Most physicists were certainly not aware of 
it and those who were aware probably doubted that infinite average values {x) or second 
moments (x 2 ) could make sense in a real phenomenon. Note that few cases which come 
under the generalized CLT were known, but they remained isolated cases (such as the 
density of first return times r in one dimension which decays as 1/r 3 / 2 for large r's). 

In the recent years, it has been more and more recognized that the generalized CLT 
could shed an interesting light on many physical processes: random walks in solutions 
of micelles [ QBL9C ], turbulent and chaotic transport [ 5ZK93| , SWS93|] , glassy dynam- 
ics [ 5hl88 , BoD95], diffusion of spectral lines in disordered solids [ ZuK94 |, thermodynam- 
ics [|TLS95| , [Tsa95| , |ZaA95| ] granu l ar flows flBoC97H , laser trapped ions [|MEZ96| , |KSW97 | 
... For reviews, see [ghl88j , poG90| , |KSZ96| , [Tsa97[| . 

The interest of the physics community for the generalized CLT seems to be stimulated 
by two important arguments. 

First, the phenomena obeying only the generalized form of the CLT, i.e. those with 
asymptotic power laws for f(x) with a < 2, exhibit a statistical behaviour which is 

1 For simplicity, any sum of independent events is called a 'Levy sum' below. 
If (x ) is finite, we say that f(x) is a 'narrow' probability density. 
3 We then say that f(x) is a 'broad' probability density. 



markedly different from the behaviour of the phenomena obeying the usual CLT. It is 
thus important to identify whether a physical process comes under the generalized form 
of the CLT (in the first place to avoid the use of natural but irrelevant concepts, such as 
the average value, derived from the usual CLT). This is illustrated in section 2 with the 
simple model of the Arrhenius cascade. 

Second, the generalized CLT provides an efficient tool for the quantitative study of 
some physical problems. This is illustrated in section 3 with a Levy flight theory of 
laser cooling of atomic gases. In this case, it is worth noting that the ideas derived from 
the generalized CLT have had practical consequences, leading to more efficient cooling 
strategies and to record low temperatures. 

Finally, the generalized CLT also provides a useful qualitative insight for some random 
walks even when it is not strictly valid. This is discussed briefly in section 4. 

Note that this contribution presents the point of view of a physicist and, as such, might not be 
rigorous on all mathematical aspects. 



2 The Arrhenius cascade 

We present the model of the Arrhenius cascade in section 2.1, which is shown to exhibit 
an unexpected statistical behaviour in section 2.2, analyzed with the generalized CLT in 
section 2.3. This toy model presents generic effects of the generalized CLT in physics. 



2.1 The model 



We consider a physical system placed in the potential landscape schematized in figure 1 
and submitted to thermal fluctuations. This model is called here the Arrhenius cascade. 
It is inspired from studies of disordered systems, like glasses relaxing towards low energy 
states [|Shl88|, BoD95], which obey similar equations. 




Figure 1. The Arrhenius cascade. 

The potential U (x) is a tilted random 'washboard'. It presents n local minima or 'wells', 
labelled by i, separated from the next minimum i + 1 by a potential 'hill' of random height 
Ui (the UiS are independent variables). We assume an exponential probability density 




f(U) for the hills U t : 

W) = i- e- u ' u ° , (3) 

where U$ is the average height of the potential hills. The variable x may be a position 
coordinate or any coordinate of the system. 

At any time, the physical system is trapped in one of the local potential wells. Due to 
thermal fluctuations, the system performs sudden jumps from one well to the other one 
downwards. The global tilt of the potential hill U(x) is large enough to prevent the system 
from performing upward jumps. Therefore, the system can only cascade downwards. The 
trapping time T{ in the well i is related to U{ by the Arrhenius lawQ: 

Tl = r e u ^ T ) , (4) 

where tq is a characteristic time, k-g is the Boltzmann constant and T is the temperature. 

We consider a gedanken experiment in which the experimentalist wants to know the 
average trapping time in a single well, but can only measure the time t n needed to go 
through all the n wells of the system: 

n 

in = ' ( 5 ) 

i=l 

To reduce the measurement uncertainty, he can repeat m times his measurement of t n . We 
assume that, between each measurement, the realization of the disorder changes, i.e. the 
numbers Ui, ...U n change (and thus t%, ...r n change accordingly). His estimation r exp (iV) 
of the average trapping time will therefore be: 

m y N 

T** P (N) = -'£t n = ( 6 ) 
i=\ i=i 

where N = m x n is the total number of explored wells and the Tj's are independent 
random variables defined by Eq. (^) and (Eh. 

2.2 Behaviour of the Arrhenius cascade 

Simulated measurements of T exp (iV) are represented in figure 2 for two different temper- 
atures. For a temperature T = 3C/o/A;b, T exp (iV) converges nicely to the average value 
(x) when iV increases, after exhibiting reasonable fluctuations at small N. This is the 
expected, standard behaviour. 

For a lower temperature T = 0.5C/o//cb, the behaviour of r exp (N) is markedly differ- 
ent. The measured T exp (Nys do not seem to converge towards any constant value but 
rather to diverge in a very fluctuating way with increasing N. A detailed analysis would 
reveal that the t„'s also fluctuate very much from one measurement to the other. This 
unusual statistical behaviour would puzzle most experimentalists: large fluctuations and 
irreproducibility are usually considered as the indication of a problem in the experimental 
procedure, arising from poorly controlled parameters. But here, one would find that the 
experimental setup works apparently well 

4 In a realistic model, tj is not deterministically fixed by Ui and the expression (^]) gives only the average 
value of Tj. Taking into account the fluctuations of n for a given Ui would not change qualitatively the 
discussion presented here. 

5 A closely related situation has recently appeared in a quantum tunneling problem. See section 4. 
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Figure 2. Simulated measurements of the estimated T cxp (N). We have chosen 
n = 10 and tq = I. The quantity r oxp (iV) for the case k^T = 3Uq has been 
multiplied by 30 to be more visible. 



2.3 Application of the generalized CLT 

Let us see how the generalized CLT can shed light on the previous observations. We 
first easily calculate the probability density f(r) of the trapping times. It is given by the 
relation /(r)dr = f(U)dU, which leads immediately to 



f(r) 



a 



-l+a 



with a 



k B T 
U 



(7) 



Having a power law for /(r), the probability densities of r exp (A r ) at large ./V are provided 
for any a by the generalized CLT. 

In the high temperature case (T = 3Uo/k^), we have a = 3 which is larger than 2. 
Thus (t 2 ) is finite and the usual CLT applies: T exp (N) is gaussianly distributed at large iV 
and tends to (r). This agrees with the observations (figure 2) and does not need further 
discussion. 

In the low temperature case (T = 0.5Uo/k B ), on the contrary, we have a = 0.5 which 
is smaller than 1. The usual CLT does not hold anymore and the specific features of the 
generalized CLT will play a crucial role. The generalized CLT tells us that we should 
consider the quantity Zn = Yli=i Ti/N 1 ^ and that the density /(Zjv) tends to a Levy law 
L q (Zjv) at large N. This theorem has several important physical consequences: scaling 
of the Levy sums, domination of the Levy sums by a single term and large fluctuations of 
the Levy sums (most of these consequences are presented in [ BoG90| ] ) . We only treat here 
the case a < 1 (and not 1 < a < 2) for which the consequences of the generalized CLT 
depart most strongly from the ones of the usual CLT. These consequences are: 



• The most probable value of the sum T/v = J2iLi T i scales as 

T N ~ N l ' a , 



(8) 



and not as N which is more usual. Practically, this implies that the time t n spent in 
an Arrhenius cascade of n wells does not scale with the size n of the cascade, but more 



rapidly due to the generalized CLT^. Similarly, the experimental value T exp (N) = 
— h~ y]TLi t, n does not tend to a constant for a large number of measurements m 
but diverges with m, as m~ l+l l a . This explains the observed behaviour in figure 2. 
Thus, the size of the system and the number of measurements play a non trivial role 
in the measured values, an unusual situation in physics. 

The notion of average value is irrelevant here since (r) = 00. One can somehow 
replace it by the notion of typical value, i.e. most probable value. When a < 1, 
the typical terms Xi of a Levy sum = x i are n °t ah of the same order (as 
they are when f(x) is a narrow density) but present a hierarchical structure. In 
particular, the typical largest term x max of the sum can be shown to be of the same 
order of magnitude as the sum itself 

N 

X N = Xi - Xmax ' ( 9 ) 
1=1 

however large N might be (Eq. (^) is valid within pref actors that do not depend on 
N ). This domination of a sum by a single term, or by a small number of terms, is a 
signature of Levy statistics in a physical problem (see figure 3 below). 

It can be used cleverly in physical experiments as a 'Levy microscope': by mea- 
suring a macroscopic quantity Xn = Ylf=\ x ii one ma y obtain an easy access to a 
microscopic information, x mSuX , while the direct measurement of x max (i.e. in the 
Arrhenius cascade, the direct study of a single well) can be physically impossible. 
Such advantageous use of the statistical domination of a single term has been made 



implicitly, for instance in studies of quantum tunneling [ R,oB84 | 



A Levy sum Xn = Yli=i x i fluctuates as much as a single term, when a < 1. This 
is a direct consequence of the domination of the sum by a single term. It can also 
be seen as a consequence of the fact that the tail of the Levy law L a (Zff), which 
determines the fluctuations, decays as l/Z^J a , exactly as the tail of f(x). This 
explains the highly fluctuating r exp (N) obtained in figure 2 as being intrinsically due 
to the type of involved statistics and not to some technical experimental problem. 
Thus, fluctuations do not vanish as usual with the increase of the size N of the 
statistical sample. The sums Xn retain an intrinsically large irreproducibility. This 
is in contradiction with a traditional motivation for applying statistical methods: to 
go beyond the irreproducibility of individual events in order to obtain quasi-perfect 
reproducibility for large ensembles of events. However, the generalized CLT still 
allows for some predictibility in the statistical sense, since it predicts the stable form 
of the probability density of Z^. 

Physics has incorporated two new types of randomness during this century: quan- 
tum uncertainty and deterministic chaos. It seems to us that the non-averaging 
out of fluctuations in Levy flights can also be recognized as an important type of 
randomness^]. 



Such 'anomalous' scaling with the system size appears in some complex phenomena like phase transi- 
tions near a critical point. What is striking here is to obtain such scaling in a very simple problem. 

7 J.P. Bouchaud speaks of a 'science of irreproducible results'. B. Mandelbrot uses the term 'wild 
randomness'. 



3 Subrecoil laser cooling of atomic gases 



In the Arrhenius cascade (section 2), a sum of iV independent terms was directly measured 
and the generalized CLT could be applied directly to analyze the results. In this section, 
we proceed a step further by studying a richer physical problem, called subrecoil laser 
cooling. In this case, Levy sums — and their properties dictated by the generalized CLT — 
play an essential role, although they are not measured directly. 

We introduce subrecoil laser cooling in section 3.1. In section 3.2, we show that power 
law densities of time variables appear, which implies the non-ergodicity of the process. 
In section 3.3, the generalized CLT is used to get some insight on the cooling efficiency. 
In section 3.4, quantitative predictions are derived, using the 'sprinkling density'. In sec- 
tion 3.5, we indicate how the insight provided by the generalized CLT enables to optimize 
the cooling strategy. 

The starting point of the approach presented here has been presented in [BBE94] and 
| Bar9E 1 ■ A detailed description of the theory will appear in |BB A99|1 . 



3.1 Subrecoil laser cooling 

Laser cooling of atomic gases is based on the momentum exchanges between photons and 
atoms. In standard (not subrecoil) laser cooling, laser configurations and atomic transitions 
are carefully chosen so that these momentum exchanges lead to a friction force. This 
friction force damps the thermal atomic momenta p, thereby reducing the momentum 
spread (standard deviation) 5p of the atomic gas, which is equivalent to reducing the 
effective temperature T defined by 

k B T = 5p 2 /M (10) 

where k-g is the Boltzmann constant and M is the mass of the atoms. Temperatures 
commonly achieved in the last ten years are in the range of a few microkelvins, 8 orders of 
magnitude below room temperature. This has opened exciting new possibilities for atomic 
and quantum physics and has been a key ingredient in the realization, in 1995, of a new 
state of matter called Bose-Einstein condensate. The 1997 Nobel prize of physics was 
attributed to S. Chu, C. Cohen- Tannoudji and W. Phillips for their contributions to laser 
cooling. 

Standard laser cooling mechanisms are fundamentally limited to temperatures larger 
than the so-called 'recoil temperature'. Indeed, among the momentum exchanges between 
atoms and photons, some — the ones due to spontaneous emission — occur in a random 
direction. Each spontaneous emission of a photon by an atom thus results in an uncon- 
trollable random recoil of the atomic momentum p by a quantity Kk, where Uk is the 
momentum of a single photon. Therefore, the standard deviation bp of the atoms is ex- 
pected to be always larger than Tik. This implies (cf. Eq. (|i~0D) laser cooling temperatures 
T larger than the recoil temperature defined by Tr = (fik) 2 / (k^M). The recoil tempera- 
ture Tr is on the order of one microkelvin for the configurations frequently used. 

To sum up, the randomness of spontaneous emission, which is essential for the cooling 
since it provides a dissipative contribution to the atomic evolution, is also harmful to 
the cooling since it implies a limit temperature. As spontaneous emission of photons by 
atoms placed in laser light seems unavoidable, the recoil temperature was for some time 
considered as an absolute limit for laser cooling. 



Subrecoil laser cooling, i.e. T < Tr, is however possible. Indeed, although spontaneous 
emission of photons is an intrinsically random quantum process, it can be partly controlled. 
The key idea is to create a spontaneous emission rate R(p) which depends on the atomic 
momentum p (figure 3) and which vanishes at p = 0. This was first proposed and realised 
in 1988 [ AAK88; , [AAK89| ] using a nice quantum effect called a 'dark resonance', because 
a resonance occurs at p = which prevents the spontaneous emission. Today, record 
low temperatures reached experimentally with dark resonances approach Tr/1000, which 
corresponds to a few nanokelvins only. 




-bp bp p 



Figure 3. Principle of subrecoil cooling, a) The spontaneous emission rate 
R(p) vanishes at momentum p — 0. b) The atoms perform a momentum 
random walk and accumulate in the vicinity of p = 0. 



One can follow the evolution of an atom with an initially non-zero momentum p. The 
spontaneous emission rateQ R(p) being large, the atom will spontaneously emit photons^] 
and therefore its momentum will change in a random way. This random walk will eventu- 
ally lead by chance the atom in the vicinity of p = where the spontaneous emission rate 
is very small (see figure 3). There, the atom stops exchanging momentum with photons 
and it remains so-to-speak 'trapped' in what is called a 'dark state'. The time r(p) of 
residence at momentum p is 

= W) • 

the time interval between two spontaneous emissions. If this residence time is long enough, 
the atom keeps the same small momentum till the end of the experiment. If not, it emits 
a spontaneous photon, which starts a new momentum diffusion process and gives to the 
atom another chance to reach the vicinity of p = 0. Thus atoms accumulate in the vicinity 



The spontaneous emission rate R(p) can be simply seen as a diffusion coefficient that has the peculiarity 
of varying with the momentum p. 

9 Bcfore each spontaneous emission, the atom absorbs a laser photon. The recoil effects of photon 
absorption, which is a deterministic process, are not essential here and are therefore ignored. 



of p = in long-lived states: a cooling effect occurs. This cooling relies on a random walk 
of the atomic momentum, unlike standard laser cooling which rests on friction forces. 



The most important question for a cooling process is to determine the typical momen- 
tum 5p at the end of the random walk or, equivalently the temperature T (cf. Eq. (|n])). It 
is difficult to answer it with the usual analytical or numerical methods of atomic physics^ 
because very different momentum and time scales are present in the problem. A conjec- 
ture was proposed in 1988 in which the interaction time 9, defined as the time the atoms 
interact with the lasers, plays a crucial role. Consider an atom reaching a momentum p 
such that the residence time r(p) is larger than the interaction time 0. This atom will 
thus keep its momentum p till the end of the experiment and will be detected with this 
momentum p. The conjecture consists in assuming that only the atoms such that 

rip) > e (12) 

keep their momentum till the end of the experiment. Obviously, this can not be strictly 
true: some atoms will reach a small momentum p after a significant time t has elapsed 
from the beginning of the interaction with the lasers so that, for them, the condition to 
stay at momentum p would rather be r(p) > 9 — t. But let us assume that condition ( |l2|) 
is the relevant criterion for the trapping of atoms. This predicts that a momentum peak 
will form with a width 5p$ given by 

T(5p e ) ~ 9 ■ (13) 

Moreover, it can be shown that the residence time r(p) varies as 

rip) cxl/p 2 ■ (14) 

Introducing this relation into Eq. ( |T3| ) gives the conjectured momentum scale 5pg which 
is reached after an interaction time 0: 

&Pe oc (15) 
or, for the temperature Tq oc (5po) 2 (cf. Eq. (|T0|)), 



Te oc i • (16) 

This result is both interesting and surprising. Interestingly, it predicts that the tem- 
perature can be reduced towards lower and lower values when the interaction time 9 is 
increased. The recoil temperature limit, which arises in standard laser cooling from spon- 
taneous emission, is no more a limit here. Here indeed, spontaneous emission is present 
to create a random walk that brings the atoms to p ~ 0, but spontaneous emission stops 
when the atoms reach a small enough momentum. Surprisingly, there does not seem to be 
any limit for the cooling. 



10 Ho wever, in one particular case, analytical solutions based on th e usual methods have been 
found |A1K96, 3SY97|. Numerical approaches have also been developped | |CBA91 ] using a new type of 



quantum simulations. 



This has motivated a series of experiments with longer and longer interaction times 
9. The recoil limit was first overcome in 1988, reaching Tg ~ Tr/2 [|AAK88j 1 . Longer 
interaction times allowed to reach Tg ~ T R /40 in 1994 HBSL94| , |Bar9E| 1 and T e ~ T R /800 



in 1997 [SHK97], establishing a new temperature record each time. Moreover, these ex- 
periments agree well with the conjectured temperature dependence of Eq. (Jl6|) . 

However, one obviously needs a better understanding of what determines the temper- 
ature Tg. A related question is the proportion of cooled atoms: can a random walk with 
no driving force lead to an accumulation of all the atoms in the vicinity of p = 0? Is 
there rather only a small proportion of cooled atoms? How does this fraction vary with 
the interaction time 9 and with the number of dimensions of the random walk? As we 
will see below, random walk techniques and the generalized CLT provide answers to these 
questions. 



3.2 Trapping time densities and non-ergodicity 

Recently developped quantum simulations [pCM92| , |DZR92j , |Car93| , |CBA91| allow to 
follow the momentum random walk of a single atom in the process of subrecoil cool- 
ing [ CBA91 , BBE94 , Bar95|| . An example of such a random walk is represented in fig- 
ure 4. We see how the random evolution of the atomic momentum sometimes leads to 
p ~ states where the atom remains for a long time because the spontaneous emission 
rate R{p) = l/r(p) vanishes in p = 0: this is the principle of subrecoil cooling at work. 
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Figure 4. (a) Example of a momentum random walk resulting from a Monte- 
Carlo simulation of subrecoil cooling of metastable helium atoms. The unit of 
atomic momentum p is the momentum Tik of the photons. The zoom (b) of 
the beginning of the time evolution is statistically analogous to the evolution 
at large scale, a fractal property typical of a Levy flight. 



More importantly, figure 4 presents an interesting statistical feature that triggered the 
Levy flight approach of laser cooling. The single largest residence time r max amounts to 
70 % of the total time 9 while the atom has occupied 4000 different momentum states 



during this total time. Thus a single event dominates a sum of a large number of events, 
which is an indication of a Levy flight. 

Is there really a Levy flight in the problem? Let us estimate the probability density f(r) 
of trapping times r, defined as the times spent by an atom in an interval [— ptrap! +Ptrap] 
called the " trap" of size ptrap smaller than the size frk of a random momentum stepp 1 ]. The 
probability density p(p) for an atom reaching the trap to fall in a state of momentum p 
can then be considered as independent of p (in one dimension) : 

- ' (17) 

The trapping time density /(r), given by /(r)dr = p(p)dp, is easily obtained using 



Eq. (|TJ): 



/(r) oc -1- with a= X -- (18) 

Thus, if we consider as a first step that the interaction time 9 is the sum of trapping 
times Tj with the density ( |l8l)Pl , the time 9 is indeed a Levy sum, whose behaviour is 
dictated by the generalized CLT with a < 1. We have a Levy flight in time, which imme- 
diately accounts for the domination of a single trapping event in figure 4 (see section 2.3). 

There is a deep physical consequence of this Levy flight, the absence of ergodicity. The 
ergodic hypothesis, an important ingredient in statistical physics, is the assumption that 
time averaging of a physical quantity yields the same result as ensemble averaging. Time 
averaging requires following a particle over a time much larger than all characteristic times 
of the problem. This is impossible here. Indeed, as the time 9 gets larger, larger trapping 
time scales (up to 9) appear and the time averaging procedure does not converge. This is 
reflected in the fact that we have a Levy flight on a time variable 7q, with infinite average 
trapping times. Thus, subrecoil cooling is a non-ergodic process. The non-ergodicity is 
associated to the absence of cooling limits. The cooling goes on for ever because larger and 
larger trapping times r(p), corresponding to lower and lower momenta p, can be reached 
with increasing 9. 



3.3 Trapping, recycling and the generalized CLT 

The history of an atom over a time 9n can be seen a series of N trapping times Tj 
interrupted by N times fj spent out of the trapP|. The times fj are the usual 'first return 
times'. We also call them 'recycling times' because the atoms coming out of the trap are 
given another opportunity to reach the trap, they are 'recycled'. 
Thus, the interaction time 9n writes as 

9 N = Tl+h + - + T N + T N = T N + f N , (19) 

11 Under these conditions, the trapping times r(p) are simply the residence times r(p) of Eq. (jTl|) and 
Eq. (jl|) in the region [— ptrapj +Pt*ap]. 

12 We neglect here the times f% spent outside the trap. These will be taken into account in the following 
section. 

13 Thus, the same non-ergodic properties occur for the Arrhenius cascade at law temperatures (see sec- 
tion 2). 

14 Note that the initial problem is a momentum random walk, which we treat efficiently by considering 
the associated random walk in time, a standard method for these problems. 



where 



is the total trapping time, and 



JV 

p n = ]T n (20) 

i=l 



N 

Tv = ]T f* (21) 

i=l 

is the total recycling time. Both T/v and T/v are sums of independent variables. The 
application of the generalized CLT to these sums gives in a simple way a qualitative 
answer for the proportion of cooled atoms, as we discuss now. 

Consider first the case in which the spontaneous emission rate R{p) tends to a non- 
vanishing constant at large p. Then, at large p, we have a standard random walk with 
a constant diffusion rate. For a ID problem, the probability density /(f) of first return 
times f is known to decay at large f as 



/(f) oc -r^-r with a = - • (22) 



It thus decays exactly in the same way as /(r). According to the generalized CLT, for 
large N's, the sums T/v and T/v behave as 

T N ~ N 1/a = N 2 , (23) 
f N ~ N 1/& = N 2 ■ (24) 

Therefore, for long times (cf. large iV's), one has T/v ~ T^: the atoms spend a finite 
fraction of their time in the trap and a finite fraction outside the trap. We thus expect the 
proportion of cooled atoms to tend to a constant, strictly between and 1. More elaborate 
calculations confirm this non-trivial result. 

Consider now the case in which a friction mechanism is added to prevent the atoms to 
diffuse to too large momenta p. This friction confines the momentum diffusion in a finite 
zone. In this case, /(f) is a narrow probability density with a finite average value. Thus, 
according to the usual CLT: 

Tn ^ N{f) ■ (25) 



Comparing this to Eq. (23) which is still valid here, one has 

T N » f N (26) 

for large N . This implies that all the atoms will be cooled, which is again confirmed by 
more elaborate calculations. 

More complicated cases can be considered by including diffusion in 2 or 3 dimensions 
or by including the 'Doppler effect' which modifies the rate R{p) at large p. In each case, 
the generalized CLT provides the asymptotic behaviours of the sums T/v and T/v which 
yield immediately the qualitative asymptotic proportion of cooled atoms. 



3.4 Momentum distribution 



Up to now, we have presented mostly qualitative results. We want to sketch here how 
some quantitative results are obtained. 

The main features of the cooling process are given by the momentum distribution 
(probability density) V(p,0) of trapped atoms at time 9. This momentum distribution 
V(p,0) writes as an integral over the times ti at which the atoms enter the trap for the 
last time: 

V( P , 9) = p(p) I dttSRtoMe - U\p) • (27) 

J 

The quantity p{p) is the probability density for an atom entering the trap to reach the 
momentum p (in one dimension, we have seen in Eq. ( |l7| ) that p{p) = l/(2p t rap))- The 
quantity Sr,(£) called here the 'sprinkling density' is the probability density for an atom to 
return into the trap at time t, regardless of the number of times the atom has entered the 
trap before. The quantity iP{t\p) = /(r')dr' is the probability that an atom remains 
in the trap during a time longer than t (where f(r) is the trapping time density defined 
in section 3.2). 

The momentum distribution can be calculated explicitly. For instance, in a simple ID 
model with infinite (r) and finite (f), one obtains 

V(p,9) = h{6)g{ApVe) (28) 

where h(9) oc y/6 is the height of the cooled peak at p = 0, A is a constant. The function 
G(q) is given by G(q) = 1 for q < 1 and by G(q) = 1 — (1 — q^ 2 ) 1 ^ for q > 1. The width 
bpQ of V(p, 9) decays as 1/VO, which proves the 1988 conjecture of a temperature decrease 
without any fundamental limit (cf. Eq. (|l6|)). This calculation can also be done for more 
complicated cases in any dimension where it is very useful. For instance, one can study 
the influence of the exponent f3 in the spontaneous emission rate R(p) oc pr , as described 
in the next section. 

The key point to obtain the momentum distribution is the calculation of the sprin- 
kling density Sn(t). The sprinkling density is obtained relatively easily with a Laplace 
transform^ The result is interesting. If (r) and (f) are finite, then S^it) tends to a 
constant at large times. This is an expected 'ergodic' result: the rate of return events is 
asymptotically constant]^]. On the contrary, if (r) or (f) is infinite, then Sn(t) decays to 
at large times. This is a signature of non-ergodicity: at large times, the density of return 
events go to because the longer and longer Tj's or fj's which tend to appear slow down 
the diffusion. Such a process has a 'history': the measurement of Sji(t) at any time tells 
when the diffusion has started. 



3.5 Optimizing laser cooling with the generalized CLT 

A remarkable outcome of the usual CLT is that the statistical behaviour of Levy sums X/v 
at large N is determined only by two parameters, (x) and (x 2 ). Thus, the detailed features 
of f(x) can be forgotten if one is interested only in the large N properties of the Levy 
sums. Similarly, with the generalized CLT in the cases a < 2, only the asymptotic power 

15 In fact, the generalized CLT is not explicitly used in the derivation of Eq. (p^). 
16 In a Poisson process, this rate is constant at any time scale. 



law behaviour of f{x) is relevant to determine the behaviour of Levy sums at large N. For 
a positive variable x, for instance, this power law is described by two parameters only, the 
exponent a and the prefactor of the power law. This can provide a useful insight when 
confronted to a complex physical phenomenon with many parameters: the generalized 
CLT shows that the many physical parameters combine into only two relevant statistical 
parameters. 

Such an insight has been applied in practice to improve a subrecoil laser cooling mech- 
anism called Raman cooling KaC92 |. Raman cooling, like the dark resonance cooling 
described in section 3.1, rests on a p-dependent spontaneous emission rate R{p) analogous 
to the one in figure 3a. The main difference is that the rate R{p) results from the su- 
perposition of pulses that can be chosen nearly arbitrarily. This gives flexibility to this 
mechanism and makes it a good case study for cooling optimization. On the other hand, 
the large number of parameters (~ 30 for the initially used sequence of pulses) to be 
optimized makes it necessary to find simplifying guidelines. 

By carefully using the generalized CLT, we have proposed a new very simple sequence 
of pulses [FtBB95|: it relies on 4 pulses only (compared to 14 initially); the shape of the 
pulses is the simplest possible while the initially used pulses were sophisticated Blackman 
pulses. 

The results are eloquent. With the initially used sequence of 14 pulses, the temperature 
Tq varied as Tq oc \j\[9 with the interaction time 9. With the new sequence of 4 pulses, 
the new shape (which changes the exponent (3 of the rate R(p) — pr from f3 = 4 to (3 = 2) 
leads to Tq oc 1/9, a much faster cooling. Moreover, with this new shape, if the pulses 
parameters (width and position) are adapted to the considered interaction time 9, one 
obtains an even faster cooling T e oc l/6> 4 / 3 [|RBB95| , |Rel96H . 

These predictions have been successfully tested experimentally and led to record low 
temperatures (2.8 ± 0.5 nK) for a cesium gas. This shows how the generalized CLT can 
have significant practical consequences. 



4 Imperfect Levy flights 

We have presented in sections 2 and 3 two examples where the generalized CLT applied 
perfectly. However, there are many physical cases where the generalized CLT is useful 
although the conditions to apply it are not, strictly speaking, mathematically fulfilled. 
This may occur either because the asymptotic decay of f(x) is not purely a power law or 
because f(x) is a truncated power law. 

Let us first discuss the truncation problem^, i.e. the cases in which f(x) decays as 
l/x 1+a for x < xq and decays more rapidly for x > xo so that (x) and (x 2 ) are finite. In 
the mathematical sense, the usual CLT applies. However, due to the power law tail, the 
convergence to the asymptotic gaussian for the probability density of the Levy sums can be 
extremely slow, being reached for N typically of 10 3 or larger [ MaS94j| , while in most cases 



for which the usual CLT applies, the approximate convergence to a gaussian is obtained 



17 In section 3, the sums TV = 53 =i Ti were limited by the available interaction time 8 which is also 
a kind of truncation. However, this truncation of the sum itself by an experimental parameter (here the 
interaction time, in other cases the system size) does not prevent the appearance of all the important 
effects of the generalized CLT; on the contrary, the fact that the truncation value, however large it may 
be, has an effect on the measured value is a signature of the generalized CLT. The truncations dealt with 
in section 4 bear on the density f(x) itself and imply a departure from the generalized CLT. 



very rapidly, with typically N ~ 4 — 5. As, in practice, one often deals with sums of a 
moderately large number N of terms, the behaviour of the Levy sums is often dictated 
by the Levy laws for relevant N values, while the gaussian behaviour is recovered only for 
irrelevantly large N values. 

Second, there are broad probability densities which decay only approximately as power 
laws. An example is provided by broad lognormal distributions, which have of course a 
finite second moment. They can be rewritten as power laws l/x 1+a ^ x ' ) with a logarith- 
mically varying exponent a(x). If the logarithmic part of a(x) is small enough, then the 
generalized CLT gives at least some qualitative guidelines for the behaviour of the Levy 
sums. We have used such guidelines to study the tunneling of electrons through a thin 
layer of insulator, a problem which has both basic and applied interests. The striking 
finding related to the generalized CLT has been that the typical current density varies by 
more than 200 depending on the scale at which it is measured [ Bar97| , DBB98| , DHB98| 
(see also [ LaB93| j), while the typical current density should be scale independent if there 
were no tails in the probability density of the current. 
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